Bloch-Zener oscillations across a merging transition of Dirac points 
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Bloch oscillations are a powerful tool to investigate spectra with Dirac points. By varying band 
parameters, Dirac points can be manipulated and merged at a topological transition towards a 
gapped phase. Under a constant force, a Fermi sea initially in the lower band performs Bloch 
oscillations and may Zener tunnel to the upper band mostly at the location of the Dirac points. 
The tunneling probability is computed from the low energy universal Hamiltonian describing the 
vicinity of the merging. The agreement with a recent experiment on cold atoms in an optical lattice 
is very good. 
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Introduction.- Dirac points in energy bands occur in 
special two-dimensional (2D) condensed matter systems 
[1], such as graphene [2], nodal points in <i-wave super- 
conductors and surface states of three-dimensional (3D) 
topological insulators [3]. They are fascinating instances 
of ultra-relativistic behavior emerging as low-energy ef- 
fective description of electrons in solids. Dirac points 
are band touching points that carry a topological charge, 
namely a Berry phase ±7r. In most systems, Dirac points 
occur in dipolar pairs (the so-called fermion doubling). 
Under variation of external parameters, it is possible to 
move these Dirac points and even make them merge. 
This merging signals a topological (Lifshitz) transition 
between a gapless phase with a disconnected Fermi sur- 
face to a gapped phase [4-7]. For example, a uniax- 
ial stress in graphene leads to a motion of the Dirac 
points but the merging transition is not reachable [8]. 
The quasi-2D organic conductor a- (BEDT-TTF) 2I3 is a 
good candidate to observe this transition under pressure 
[9], however it has not been realized yet. 

Recently, a new type of experiment, realized with ul- 
tracold atoms loaded in a 2D optical lattice, has provided 
an alternative way to study Dirac points [10]. By combin- 
ing techniques of Bloch oscillations and adiabatic map- 
ping of cold atoms, the band structure of the system can 
be studied with momentum resolution [11-14]. Specif- 
ically, the experiment of ETH Zurich [10] utilizes such 
techniques for a non-interacting Fermi gas in a tunable 
two-band system featuring Dirac points. Their existence 
is revealed through Landau-Zener (LZ) tunneling from 
the lower to the upper band. As the lattice amplitude is 
varied, a drastic change in the transferred atomic frac- 
tion provides a qualitative signature of the Dirac points 
and their merging. 

In this letter, we present a complete description of 
Landau-Zener tunneling through a pair of Dirac points, 
using a universal low-energy Hamiltonian describing the 
merging transition [7]. We show how the transferred 
fraction provides a key signature of the merging tran- 
sition and that it depends crucially on the direction of 
the motion with respect to the merging direction. We 
find a very good agreement between the computed aver- 




FIG. 1: (color online), (a) Square lattice indicating the hop- 
ping amplitudes and the two inequivalent sites, (b) Band 
structure in the gapless D phase for t' — t — 0.2, t" — 0.05 in 
units of Er, see text. The first Brillouin zone is indicated by 
the square. 

aged LZ probabilities and the experimental data. Fur- 
thermore, new experimental signatures for varying Bloch 
oscillations and a coherent Stiickelberg interferometry are 
presented. 

Tight-binding model- We consider a nearest-neighbor 
tight-binding model on a square lattice. The four hop- 
ping amplitudes between neighbors are taken as t, t along 
y- and t' : t" along x-direction (see Fig. 1). When t' ^ t" 
there are two inequivalent sites - called A and B - per 
unit cell giving rise to two bands. When t" = 0, a link is 
broken realizing a brick- wall lattice, which has the same 
connectivity as the honeycomb lattice albeit with a rect- 
angular geometry. When t' = t" , it is a standard square 
lattice with anisotropic amplitudes along x and y and a 
single site per unit cell. The Hamiltonian (with nearest 
neighbor distance a = 1 and Ti = 1) reads 

H = ( /*(k) /( o k) ) W 

with /(k) = -(te iky + te~ iky + t' e ikx + t"e~ ikx ) where 
the hopping amplitudes are positive and k = (k x ,k y ) 
is the Bloch wavevector. The energy spectrum is given 
by e(k) = ±|/(k)|. It features two Dirac cones when 
t' + t" < 2t (gapless D phase) and a gap when t' + t" > 2t 



(gapped G phase). At t' + t" = 2t the Dirac points merge 
at k = (0, 7r) and there is a single touching point between 
the two bands. When t' = t" the band structure is that 
of a square lattice, with lines of Dirac points (L phase), 
which becomes isotropic when t' =t"=t (I phase) . This 
model can therefore describe the transition between the 
G and D phase, and moreover, the crossover from the D 
to L phase. In the following, energies are measured in 
units of the recoil energy = n 2 Ti 2 / (2ma 2 ) where m is 
the atomic mass. 

Mapping to the universal Hamiltonian.- In a crystal 
which is time-reversal and inversion symmetric, merg- 
ing can only occur at G/2 points where G is a recip- 
rocal lattice vector. Near such a point, it is possible to 
write a minimal low energy Hamiltonian that captures 
the topological transition and describes both Dirac cones 
at once [7] . Close to the merging, an expansion for small 
q = k — G/2 gives rise to an effective Hamiltonian which 
has the universal form 

H= ( 2 ° A * + 2^ - \ ( 2 ) 

\A* + ^L_+ic x q x / 



and a spectrum e = ± y ( + A*) 2 + c 2 q 2 . The 
model depends on three independent parameters A*, m* 
and c x . The merging transition is driven by the pa- 
rameter A* hereafter called the merging gap. When 
A* < 0, the spectrum contains two Dirac points at 
= (0, ±^2ra*|A*|) and |A*| represents the energy 
of the saddle points connecting them, which is located 
at q = (D phase). When increasing A* towards 0, 
the two Dirac points approach each other along the q y 
direction until they merge when A* = 0. Exactly at 
the merging, there is a single touching point between 
the two bands, with a semi-Dirac (or hybrid) spectrum 

e = ±<J(q 2 /2m*) 2 + (c x q x ) 2 [5, 15]. By increasing the 

driving parameter still further, a true gap of magnitude 
2 A* > opens at q = (G phase). We map the tight- 
binding to the universal model by comparing their energy 
expansions near the Dirac points [7]. In the D phase, we 
find A* = t' + t" — 2t and the Dirac cone velocities are 
Cx — t' — t" and c y = \/4:t 2 — (t f + t") 2 , so that the mass 
is obtained from m* = -2A*/c 2 = 2/{2t + t' + t"). In 
the G phase, A* and c x are unchanged and m* = l/(2t). 

Landau-Zener tunneling with Dirac cones.- Consider 
atoms initially in the lower band performing Bloch os- 
cillations under the influence of a constant applied force 
F. By accelerating these atoms in the vicinity of a Dirac 
point, their tunneling probability to the upper band is 
finite, a problem considered by Landau and Zener [16]. 
In the following, the universal low-energy Hamiltonian is 
used to compute the interband transition probability. 

Motion along the k x direction: single Dirac cone. - In 
the D phase, an atom moving along the k x direction en- 
counters at most one Dirac cone during a single Bloch 
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FIG. 2: (color online). Motion along k x . (a) LZ probability 
Pf (eq. (3)) as a function of — A*/P and of the transverse 
momentum q y . Here c x /F « 0.5 and m*F « 0.13. (b) Tra- 
jectories along the k x direction, (c) Transferred fraction to 
the upper band t; x as a function of — A*/P for different sizes 
kpy of the cloud. 

oscillation (Fig. 2b). The LZ probability for such a lin- 
ear avoided band crossing is given by [16] 

n (gap/2) 2 t (^+A,) 2 
Pf = e C * F =e c xF (3) 

where q y is the position with respect to the merging point 
G/2 (0,tt). Note that Pf = 1 for q y ±q D = 
±\/—2m*A*, positions of the two Dirac points. Actu- 
ally Eq. (3) is not only valid in the D but extends to 
the G phase across the merging transition. This quan- 
tity is shown in Fig. 2a as a function of the transverse 
momentum q y and A*. 

As the experiment is performed with a cloud of non- 
interacting fermions, we need also to average the LZ 
probability over the initial distribution of atoms. We 
consider a 2D cloud of harmonically trapped fermions 
at zero temperature for a filling fraction sufficiently 
smaller than half-filling. The energy spectrum close 
to k = is e(k) « k 2 /(2m x ) + k 2 /(2m y ) (as mea- 
sured from the band bottom) with the band masses 
m x = (2t + 1' + t")/[U't" + 2t(t! + t")\ and m y = l/(2t). 
The semiclassical energy of an atom is therefore e(k, r) = 
k 2 x /(2m x ) + k 2 l(2m y ) + (m x uj x x 2 + m y uj 2 y 2 )/2 where 
u x /2tt and u y /2ir are the trapping frequencies [17]. The 
fraction £ x of atoms transferred to the upper band is then 
given by the averaged probability £ x = (Pf ) where 

J e{Kr)<eF dk x dk y dxdy- 

Je(k,r)<e F dk x dk y dxd V 

and 6p = k 2 Fx /(2m x ) = kp /(2m y ) is the Fermi energy, 
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which defines kp x an d ^Fy [18]. The transferred fraction 




FIG. 3: (color online). Motion along k y . (a) Total probability 
Pt for atoms tunneling to the upper band (eq. (6)) as a 
function of c x /F and of the transverse momentum q x . Here 
A*/F = -5 and m*F « 0.13. (b) Double LZ tunneling along 
the k y direction, (c) Transferred fraction £ y as a function of 
c x /F for various sizes kFx of the initial cloud. 

£ x as a function of A* and of the size kp y of the cloud is 
shown in Fig. 2c. For a cloud of finite size kp y , only a 
finite proportion of atoms may tunnel to the upper band 
when A* < 0. 

Motion along the k y direction: double Dirac cone.- 
In the G phase, the tunneling probability is vanishingly 
small. In the following, we concentrate on the D phase 
where atoms performing one Bloch oscillation in the k y 
direction have the possibility to encounter two inequiv- 
alent Dirac cones successively. The scenario is therefore 
richer than before since the tunneling process implies two 
successive Landau-Zener events (see Fig. 3b). The prob- 
ability P| associated with each LZ event is now 
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which defines the adiabaticity parameter 5. In the fol- 
lowing, we calculate the total interband probability 
associated with the two successive events, in the limit 
where they can be considered independent. Quantita- 
tively, the LZ tunneling time ~ max(V^, S)/c x q x [19] 
should be shorter than the time 2^/2m* |A* \/F it takes 
an atom to travel between the two Dirac points, i.e. not 
too close to the merging transition. 

First assuming that the two tunneling events are inco- 
herent, we combine the probabilities to find the interband 
transition probability [18, 19] 



2Pf(l 



Pi) 



(6) 



which is shown in Fig. 3a as a function of c x and the 
transverse momentum q x . Notice that Pj/ vanishes when 
q x = because P| = 1. For an initial cloud of size kp x , 
the transferred fraction is ^ y — (P^) where the average 
is defined in Eq. (A. 6) [18]. The result is shown in Fig. 
3c. 

Comparison to the experiment. - The ETH experiment 
is performed on a harmonically trapped 3D Fermi gas 
loaded in a 2D optical trap [10]. Our model treats a 
trapped 2D Fermi gas in a 2D band structure. The op- 
tical lattice potential is V(x, y) = —Vx cos 2 (to + 0/2) — 
Vx cos 2 (to) — Vy cos 2 (7n/) — 2ay/VxVy cos(to) cos(7n/), 
where a = 0.9, 6 = 7r, and the laser wavelength is 2a with 
amplitudes Vy = 1.8, < V> < 6.5 and < V x < 1. 
To make a precise comparison, we perform single-particle 
numerical band structure calculation provided by the 2D 
optical potential using a truncated plane- wave expansion, 
and establish a map between the optical lattice param- 
eters and that of the universal Hamiltonian of Eq. (2) 
for A* < 0. The latter is done by fitting the parameters 
of the universal Hamiltonian to the two lowest lying en- 
ergy bands in the vicinity of the Dirac cones [18]. The 
corresponding (t, t' , ^^-tight-binding model is then ob- 
tained by the analytical mapping (see paragraph after 
Eq. (2)) - the main features that are not captured are 
the Dirac cones tilting and the particle-hole asymmetry, 
which appear not to be relevant. Qualitatively, t stays 
roughly constant in the considered range of optical lat- 
tice parameters. The hopping t' increases when Vx — Vx 
increases, whereas t" increases when Vx + Vx decreases. 
Finally we compute the transferred fractions £ x and £ y as 
a function of Vx and Vx, shown in Fig. 4. For the exper- 
imental parameters F = 0.02, kp y — tt/2 and kp x — 2 
(corresponding to ep w 0.4), we find a striking agreement 
with Figs. 4a and 4b of Ref. [10]. 

Discussion.- First consider the case of the motion 
along k x . The line of maximum transfer probability ^ x 
(red region in Fig. 4a) corresponds to a maximal LZ 
probability Pg ^ 1 for a large number of atoms. Playing 
with the averaging order gives £ x « exp[— 7r((g 2 /2ra* + 
A*) 2 )/(c x F)], in which (q 2 y ) = k 2 Fy /6 and (<$) = fc^/16. 
Its maximum occurs when A* ~ — (g 2 )/(2ra*), which ex- 
plains why it is near the merging line A* = 0, but slightly 
inside the D phase, as seen experimentally. The trans- 
ferred fraction £ x is a symmetric function of its natural 
variable A* - when doing the average properly it is ac- 
tually slightly asymmetric - and its width reduces when 
decreasing c x (by decreasing Vx). Both features are seen 
experimentally, see Fig. 4a of Ref. [10]. Obviously, the 
width and the maximum should increase when increasing 
F. 

Adding a staggered on-site potential +A/2 opens a gap 
A in the spectrum. Experimentally, this gap is controlled 
by the parameter of the potential V(x, y). It was found 
that the decay of the transferred fraction £ x as a function 
of is well fitted by a gaussian (see Fig. 2b of Ref. 
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FIG. 4: (color online). a) Transferred fraction ^ x as a 
function of the optical lattice parameters V x and Vx (here 
kFy — 7r/2,F = 0.02). Inset: lines of constant A*, b) Trans- 
ferred fraction ^ y (kFx = 2, F = 0.02). Inset: lines of constant 
Cx/cy. c) Same as (b) with F = 0.1. d) Same parameters as 
(b) taking coherence into account and leading to Stiickelberg 
oscillations, e) Phase diagram in the (—A*, c x ) plane showing 
the G, D, L, and I phases (see text). P§ = 1 along the merging 
transition (continuous line). The crossover line PJ 7 = 1/2 cor- 
responds to c x oc F 1 / 2 \A*\ 1 / 4 and is plotted for two different 
forces (dashed and dotted lines). The color code for (a)-(d) is 
such that £ m = 0.5, 0.3, 0.3, and 0.6, respectively. The black 
line in (a)-(d) indicates the merging transition A* = 0. 



[10]). Here, we prove that it is of the form <^(A) = 
^(0)exp(-7rA 2 /(4c ;E F)) and that A « 4.3(0/tt - 1) for 
the experimental parameters V x xy ~ [3-6, 0.28, 1.8]. 

Next, consider the case of the motion along k y . 
The interband transition probability (Eq. 6) is a non 
monotonous function of the LZ probability P| and it 
is maximum when P| = 1/2. This explains the exis- 
tence of the maximum of £ y well inside the D phase (red 
region in Fig. 4b). When playing with the averaging 



(? x jc y ~ F In 2/ (ir (q%)) = const. Compared to the pre- 
vious case, £ y is a much more asymmetric function of 
its natural variable X cx c x /c y : the decay at large X is 
slower than at small X (the signal extends more towards 
large c Xl which is seen experimentally). In contrast to the 
previous case, the position of the line of maxima depends 
on F, as shown in Figs. 4b and 4c, but its amplitude is 
almost independent. Its width decreases when c x — > 0. 
Furthermore, £ y vanishes as Vx — >• 0, eventually reaching 
the square lattice limit (L phase). 

Since the two LZ events along k y are expected to be co- 
herent, Stiickelberg interferences in the transferred frac- 
tion £ y should be observable. Eq. (6) should indeed 
be replaced by P£ = 4P|(1 — P|) cos 2 ((/?/2 + ipd) where 

<p = 4 f« D dq y ^(A* + g2/2m*) 2 + c 2 x ql/F is a dynami- 
cally acquired phase in between the two tunneling events 
and (fd = — 7r/4+(5(ln(5 — l) + argT(l — iS) is a phase delay 
given in terms of the gamma function and S is defined in 
Eq. (5) [18, 19]. Averaging P^ over the 2D atomic dis- 
tribution gives the transferred fraction £ y shown in Fig. 
4d. However, such interferences are not observed in the 
ETH experiment, which we attribute to averaging over 
the third spatial direction. Briefly, interference fringes 
in Fig. 4d are lines of constant A* with a fringe spac- 
ing ~ O.OAEr. For the experimentally given trapping 
frequency in the z direction, we estimate that A* varies 
along z by ~ 0.03P^r. This should be enough to wash 
out the interferences. 

Conclusion.- Landau-Zener tunneling conveniently 
probes the energy spectrum in the vicinity of Dirac 
points. Depending on the direction of the applied force, 
the atoms experience a LZ transition through a single 
or a pair of Dirac cones. We calculated the transferred 
fraction in the framework of the universal Hamiltonian 
describing the merging transition, and found a very good 
agreement with the ETH experiment. To summarize, the 
important parameters are the merging gap A* and the 
velocity c x perpendicular to the merging direction. A 
simplified phase diagram is shown in Fig. 4e. Although 
the transfer through a single Dirac cone probes the merg- 
ing transition, the transfer through a pair of cones signals 
a double Landau-Zener event inside the D phase and a 
cross-over towards the L phase. 

As perspectives, we expect Stiickelberg interferences to 
be observable in the strictly 2D regime. Furthermore, it 
should now be possible to tune the optical lattice right at 
the merging transition and to study the semi-Dirac spec- 
trum, thus opening the way to explore new phenomena. 
For example, applying an artificial U(l) gauge potential 
[20, 21] should reveal unusual Landau levels [5]. 

We acknowledge support from the Nanosim Graphene 
project under grant number ANR-09-NANO-016-01. 
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(1 - e~ x )/2 where X = nc 2 x (q 2 x ) /c y F. 
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Appendix 



NUMERICAL BAND STRUCTURE 
CALCULATION 

In the experiment [Al], a 2D tunable optical potential 
of the form 

V(x,y) = -V x cos 2 (kx + 0/2) - V x cos 2 (kx) 
—Vy cos 2 (ky) — 2ol\JVxVy cos(kx) 
x cos(ky) cos(cp) (A.l) 

is utilized. Here, V x , Vx and Vy are proportional to 
the tunable laser intensities, a = 0.9 is the visibil- 
ity of the interference pattern, and other parameters 
= 7r, ip = are fixed for the most part in the ex- 
periment, k = 2tt/X is the laser wave vector, and the 
nearest neighbor lattice distance used in the main text 
a = A/2. The Bravais lattice of the optical potential 
is given by ai = (A/2) (1,-1), a 2 = (A/2) (1,1), and 
the corresponding reciprocal vector is given as bi = 
(2tt/A) (1,-1), b 2 = (2tt/A) (1,1). The single-particle 
bandstructure is obtained by numerically solving the 2D 
Schrodinger equation 



h 2 

2mo 



A + V(r) 



V(r) = £ty(r). 



(A.2) 



By writing the wavefunction in B loch's form iftk(r) = 
exp(ik • r) J2 m nez u ( m i n ) exp(zrabi • r + mb 2 • r), the 
Schrodinger equation in terms of the Fourier components 
u(m, n) becomes 

h 2 

{ o — \(k x + 27r(m + n)/\) 2 + (k y + 2it(m - n)/X) 2 } 

ZTYIq 

+5i}u(m, n) + {#2 u(m — 1, n — 1) + £3 u{m — 1, n) 
+^4 u(m — 1, n + 1) + £3 u(m, n — 1) + 63 u(m, n + 1) 
+£4 u(m + 1, n — 1) + £3 u(m + 1, n) 
+<5 2 u(ra + l,n + 1)} = Eu(m,ri), (A.3) 

where the various constants are defined as Si = — (V x + 
V x ~ V Y )/2, S 2 = (V x - Vy)/4, ^3 = -OAbVVxV^ and 
64 = Vy/A. We truncate the Fourier coefficients u(m,n) 
at |m|, \n\ = 2 and solve the resulting 25 secular equa- 
tions in Mathematica to obtain the discrete energy bands 
Ei(\s) for each k. A convergence test has been carried out 
for the numerical bandstructure calculation up to a trun- 
cation at |m|, I n I = 3, which yields an excellent agreement 
with the |m|, \n\ = 2 case. 

A typical bandstructure for a particular set of laser 
amplitudes is shown in Fig. Al(a), in the convenient re- 
coil energy unit Er — h 2 /2mo\ 2 where mo is the atomic 
mass. From the two lowest bands, see Fig. Al(b), we 
then determine the magnitude of the gap A*, the dis- 
tance between the two Dirac points qo and the slope c x 



6 




FIG. Al: (a) Bandstructure for V x = 5E R , Vx = 0.3 E R , 
Vy — 1-SEr in the k y -direction with k x = 0. (b) The two 
lowest bands featuring two inequivalent Dirac points, (c) The 
three phases D, G and L (dashed line) realized with the opti- 
cal potential. The starred point indicates the position where 
the isotropic square lattice model (I phase) is realized. 



around the Dirac points in the x direction as a function of 
Vx and Vx, for a fixed Vy = 1.8 Er. The interpolation 
formulae in the gapless phase are obtained as follow: 

\A*(Vx,V x )\ = (0.1134 + 0.0218 Vx - 0.0205 V%) 
xln[l + (Vx -3.8234 ln[2.7560 V x + 1.6325]) 1>2 ] , 
q D (Vx,V x ) = (0.65 + 0.1 V x ) 
x (Vx - 3.8234 ln[2.7560 V x + 1.6325]) C 
c,(^,Vx) = 1.1204e- - 3023 ^ 



v 0.42+0.05 V x 



xV 



0.3791 ln[1.8380 V x -0.7251] 



X 



(A.4) 



The slope in the ^/-direction around a Dirac cone is then 
given by c y = 2\A*\/qn. Furthermore, the boundary 
between the gapped and the gapless phase is given by 



Vx = ai In [a 2 V x + a 3 ] 



(A.5) 



with ai = 3.8234, a 2 = 2.7360, a 3 = 1.6325. After trans- 
lating into the more intuitive anisotropic square lattice 
model with t,t',t" hopping amplitudes (see the para- 
graph after Eqn (2) in the main text), we then locate 
the various phases that are realized as a function of the 
laser amplitudes, see Fig. Al(c). 



AVERAGING OVER THE ATOMIC 
DISTRIBUTION 

Motion along the k x direction: single Dirac cone 

As the experiment is done with a cloud of non- 
interacting Fermi gas, we need to average the Landau- 
Zener probability over the distribution of atoms. Taking 
a 2D cloud of harmonically trapped atoms at zero tem- 
perature filled up to the Fermi energy cp (measured from 
the band bottom), the fraction of atoms transferred to 
the upper band is given by the averaged tunneling prob- 
ability (Pf ) where the average is defined as 



(...) 



e(k,r)<e F 



dk x dk y dxdy • 



e(k,r)<e F 



dk x dkydxdy 



(A.6) 



Since the low energy expansion of the spectrum e(k) as- 
sociated with the tight-binding Hamiltonian (1) is e(k) = 
U'k 2 /(2t J rt') J rtk 2 (as measured from the band bottom), 
the semiclassical energy of an atom in the 2D anisotropic 

harmonic oscillator is e(k, r) = + + l(m x uj x x 2 + 

rriyUJyy 2 ) , where m x = (2b + t')/(2tt') and m y = l/(2t) 
are the band masses and uj x /2it and uj y j2 r K the trapping 
frequencies [A2], the average is given by 



(...) 
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f 

Jo 



dk 



yK^Fy 



2\3/2 



(A.7) 



h z 1.2 

2m y 2m x 2m y 

defines the 2D Fermi surface. 



where ep 



\(m x uj x x 2 + m y uj 2 y 2 ) 



Motion along the k y direction: double Dirac cone 

For an initial cloud of size kp x , the transferred fraction 
is (Pt) where the average is defined in (A.6) and becomes 



(...) 



16 



dk x (kj 7x k 2 ) 3 / 2 . . . (Ai 



where kp x — V^ m x e F- 



STUCKELBERG OSCILLATIONS 

In the following, we calculate the total inter-band prob- 
ability associated to the two successive LZ events, in the 
limit where they can be considered to be separated. This 
is the case if the system is in the D phase and not too 
close to the merging. Quantitatively, the Zener tunneling 
time [A3] r ~ max(y/5 , 6) / c x q x should be shorter than 
the time T = 2q D /F = 2y/2m*\A*\/F it takes to travel 
between the two Dirac points, where S = c 2 q 2 /(2c y F) is 
the adiabaticity parameter. 



7 



If the sequence between the two tunneling events is 
coherent, amplitudes have to be considered instead of 
probabilities, effectively realizing a coherent Stiickelberg 
interferometer [A3]. In this case, the total probabil- 
ity amplitude to go from the lower to the upper band 
is the sum of the amplitude for two distinct path. In 
path 1, the atom jumps to the upper band at the first 
avoided band crossing (Dirac cone) and then stays in the 
upper band at the second, such that the amplitude is 
A 1 = v/Pfe^ x e^ 1 x yjl - Pf where y/^e itpd is 
the amplitude to jump at a single avoided band cross- 
ing - the associated probability of a single LZ event 
being Pf - and tpi = dtE upper hand (t) is the phase 
dynamically acquired by the atom traveling in the up- 
per band between the two Dirac cones. The phase de- 
lay (fd which is accumulated at each tunneling event is 
given by cp^ = — 7r/4 + S(\nS — 1) + argr(l — iS) where 
T(z) is the gamma function [A3]. Up to a it/ 2 shift, 
this phase is the so-called Stokes phase. In path 2, the 
atom stays in the lower band at the first Dirac cone and 
then jumps to the upper band at the second. The asso- 
ciated amplitude is A 2 = y/l - Pf x e^ 2 x y/P^e~ i(pd 

where if2 = f Q T dtE\ ower band (t) is the dynamically ac- 
quired phase of the atom traveling in the lower band 
from one Dirac cone to the other. The total probabil- 
ity is therefore 

P? = \A 1+ A 2 \ 2 = 4Pf(l - Pf) cos 2 (^/2 + <p d ) (A.9) 

where (p = (pi — (p 2 is the dynamically accumulated phase 
between the two tunneling events 

<p = / AE(t)dt=- / AE(q y )dq y (A.IO) 

JO * J-q D 



with AE = ^ U pper band - ^lower band- For the universal 
Hamiltonian, we have 



4 /V 2 ™*I A *I // Q 2 \^ 

=fL ^v(^ +a *) +clql (A - n) 



which can be written in the form 



= 4y /2^\^!l I (^^j (A.12) 



where the function I(x) is given by I(x) = 
Jq 1 duyj (u 2 — l) 2 + x 2 and well approximated by I{x) ~ 
^4/9 + x 2 . From Eq. (A.9), the incoherent limit (Eq. 
(6) in the main text) is easily recovered by averaging over 
the phase. 
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